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Abstract 

The two-dimensional one-component plasma, i.e. the system of point- 
like charged particles embedded in a homogeneous neutralizing back- 
ground, is studied on the surface of a cylinder of finite circumference, or 
equivalently in a semiperiodic strip of finite width. The model has been 
solved exactly by Choquard et al. at the free-fermion coupling T = 2: in 
the thermodynamic limit of an infinitely long strip, the particle density 
turns out to be a nonconstant periodic function in space and the system 
exhibits long-range order of the Wigner-crystal type. The aim of this pa- 
per is to describe, qualitatively as well as quantitatively, the crystalline 
state for a larger set of couplings r = 27 (7 = 1, 2 ... a positive integer) 
when the plasma is mappable onto a one-dimensional fermionic theory. 
The fermionic formalism, supplemented by some periodicity assumptions, 
reveals that the density profile results from a hierarchy of Gaussians with 
a uniform variance but with different amplitudes. The number and spa- 
tial positions of these Gaussians within an elementary cell depend on the 
particular value of 7. Analytic results are supported by the exact solution 
at 7 = 1 (r = 2) and by exact finite-size calculations at 7 = 2, 3. 

KEY WORDS: Two-dimensional jcllium; semiperiodic boundary conditions; 
translation symmetry breaking. 
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1 INTRODUCTION 



According to the laws of electrostatics, the Coulomb potential v at a spatial 
position r G R v of the ^-dimensional Euclidean space, induced by a unit charge 
at the origin 0, is defined as the solution of the Poisson equation 

Av(r) = -s v 5(r) (1.1) 

where s„ is the surface area of the unit sphere in R u . The pair interaction 
energy of particles with charges q and q' , localized at the respective positions r 
and r', is given by 

v{v,q\v',q l ) = qq l v{\v-v l \) (1.2) 
In one dimension (ID), s\ = 2 and the solution of ijl.lfl reads 

v(x) = -\x\, v = \ (1.3) 

In 2D, S2 = 2tt and the solution of Ql.lfl . subject to the boundary condition 
Vti(r) — > as |r| — > oo, reads 

u(r)=-ln^V ^ = 2 (1.4) 

where ro is a free length constant which fixes the zero point of the potential. 
The Coulomb potential defined by Eq. Ijl.lf) exhibits in the Fourier k-space the 
characteristic singular |k|~ 2 form. This maintains many generic properties, like 
the sum rules PP, of "real" 3D Coulomb systems with the interaction potential 
v(r) = l/|r|, r e R 3 . 

The present paper deals with the equilibrium properties of the classical (i.e. 
non-quantum) one-component plasma, sometimes called jellium, formulated in 
ID or quasi-lD domains. The jellium model consists of only one mobile pointlikc 
particle species of charge q embedded in a fixed background of charge — q and 
density n such that the system as a whole is neutral. 

Thermodynamics of the ID jellium has been obtained exactly a long time 
ago by Baxter [2]. It was proven subsequently that the ID jellium is never 
in a fluid state, but forms a Wigner crystal In particular, choosing the 

free (hard walls) boundary conditions and going to the infinite volume limit, 
the one-particle density becomes periodic in space with period 1/n. This long- 
range order is present for all densities n and all temperatures. Although the 
ID jellium is not in a fluid state, it behaves as a conductor in the sense that 
arbitrary boundary charges are perfectly screened by means of a global transport 
of the particle lattice in the background, with no additional polarization in the 
bulk 0. 

Translation symmetry breaking was documented also on a quasi-lD system, 
namely the 2D one-component plasma living on the surface of a cylinder of cir- 
cumference W [HJ ■ This system is exactly solvable at the dimensionless coupling 
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constant T = 2 7 . In the thermodynamic limit of an infinitely long cylinder, 
the one-particle density is given by an array of equidistant identical Gaussians 
along the cylinder's axis, with period l/(nW). 

In ID and quasi-lD Coulomb systems, the variance of the charge in an 
interval / remains uniformly bounded as |/| — > oo. The existence of periodic 
structures is related to this boundedness of the charge fluctuations 

The present work proceeds in the study of the 2D jellium on the cylinder 
surface [51 [7] • Our aim is to describe, qualitatively as well as quantitatively, 
the crystalline state for a larger set of couplings V = 2<y (7 = 1, 2, ... a positive 
integer). At these couplings, the underlying model is shown to be mappable 
onto a ID anticommuting-field theory following the method of Ref. [TJ]|, and 
its density profile is expressible in terms of the corresponding field correlators. 
The assumption of the periodicity of the particle density in the thermodynamic 
limit reveals uniquely that the density profile results from a superposition of a 
hierarchy of nonidentical Gaussians with a uniform variance but with different 
amplitudes. The number and spatial positions of these Gaussians within an 
elementary cell depend on the particular value of 7. The analytic results for the 
crystalline state are supported by the exact solution at 7 = 1 (T = 2) and by 
the exact finite-size calculations at 7 = 2,3. 

The paper is organized as follows. In Section 2, we present basic formulas 
for the one-component plasma living on the cylinder surface. Section 3 deals 
with the ID fermionic representation of the model for the special values of the 
coupling constant T = 27 (7 a positive integer). Section 4 is devoted to a 
general analysis of the density profile in the thermodynamic limit; the Gaussian 
structure of the crystalline state is revealed for any value of 7. The analytic 
results are verified in Section 5 on the exact solution of the model at 7 = 1, and 
by the exact finite-size calculations at 7 = 2 and 7 = 3. 



2 THE MODEL 

First we define the 2D one-component plasma confined to the surface of a cylin- 
der of circumference W and finite length L, in the canonical ensemble. The 
cylinder surface can be represented as a 2D semiperiodic rectangle domain A 
with r = (x,y) € A if — L/2 < x < L/2 (free or hard walls boundary condi- 
tions at x = ±L/2) and — W/2 < y < W/2 (periodic boundary conditions at 
y = ±W/2). It is sometimes useful to use the complex coordinates z — x+iy and 
z — x — iy. There are N mobile pointlike particles of charge q in A, embedded 
in a homogeneous background of charge density = —qn with 

n =m (2 - 1} 

so that the system as a whole is neutral. The interaction potential between 
two unit charges at 17 and r2 is given by the 2D Poisson equation (f 1 . 1|) with 
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the requirement of periodicity along the j/-axis with period W. Writing the 
potential as a Fourier series in y, one gets [H] 



u(ri,r 2 ) = - In 



2sinh 



7r(zi - Z 2 ) 



W 



(2.2) 



At small distances |ri — r 2 | << W, this potential behaves like the 2D Coulomb 
potential (|1.4f> with the constant ro = W/(2ir). At large distances along the 
cylinder \x\ — x 2 \ » W, this potential behaves like the ID Coulomb potential 
— (ir/W)\xi — x 2 |. We shall need to express the absolute value on the rhs of 
H2.2JI formally as the product g(z\)g{z 2 )\f {zi) — /(z 2 )|. This can be done in two 
ways: 



2 sinh 



n(zi - z 2 ) 



W 



e w 



e w " 



(2.3) 
(2.4) 



The final results cannot depend on the particular choice, and we shall adopt 
the representation (j2.3|l . For a given configuration {ri, . . . , rjy} of charges, the 
total energy of the particle-background system is given by '& 



E N ({r}) = q 2 v i Y 3 1 r fc) + ™ q 2 J2x 2 + B N 

j<k 3 

where Bn is the background-background interaction constant. 
The partition function at inverse temperature (3 is defined by 



(2.5) 



Zn = 



1 

ah 



N 

II d V 



(2.6) 



It depends on the dimensionless combination T — (3q 2 called the coupling. The 
multiplication of Zn by a constant does not effect the particle distribution 
functions, so for notational convenience we omit the interaction constant Bm in 
l|2.5|l and multiply each volume element d 2 rj in (|2.6|l by W~ 2 , to get 



i r N 

z * = / n [d'vo^i)] n i e#z 



e w 



'2=1 j<fe 

Here, w is the one-body Boltzmann factor 



w(z, z) = w(x) = exp 



Trr 



-^--(iV-l), 



(2.7) 



(2.8) 
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The particle density at point r € A is defined as 

JV 

n(r) = ( YA*-*i)) (2-9) 



where (• • •) denotes the usual canonical average. It can be obtained in a standard 
way as the functional derivative 

n(z,z)=w(z,z)- — 7 — TT-lnZ/v (2-10) 

ow{z,z) 

Due to the cylinder geometry of the system, the particle density depends only 
on the x-coordinate, n(r) = n(x), and exhibits the reflection symmetry n{x) = 
n{— x) with respect to the origin x = 0. The charge neutrality of the system is 
equivalent to the condition 

L/2 

dx[n(x)-n] = (2.11) 

-L/2 

Our task is to determine the particle density profile in the thermodynamic 
limit N, L — * oo (the circumference W of the cylinder is finite), where the 
background density n given by (|2.1() stays constant. The exact ID solution of 
the jellium at any temperature [2] and the exact 2D solution at the coupling 
r = 2 7 indicate two characteristic features of this density profile: 

• The thermodynamic limit of the density profile depends on which one 
of the two subsequences, the particle number N = even and N — odd 
integers, is chosen; we denote by n^ e '(x) and n^°\x) the corresponding 
density profiles defined in — oo < x < oo. The plots of rS e '°\x) are 
expected to be periodic with a period A, 

n {e '°\x± A) = n^°\x) (2.12) 

The two density profiles are supposed to be the two realizations of the 
same periodic function shifted to one another by a half period, 

n^\x±\/2) =n(°\x) (2.13) 

• The period A is equal to 1/n in ID, independently of the temperature. 
The exact 2D solution for the coupling T = 2 gives A — 1/ (nW) . In both 
cases the elementary cell with the size of the period in the x-direction 
contains just one particle 1 . This motivates us to suggest that the period 

A = — (2.14) 

is present for an arbitrary coupling T. 



we are grateful to B. Jancovici for noticing to us this important fact 
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These two working hypothesis will be incorporated into an analytic treatment 
of the model, and subsequently justified numerically with a high precision by 
finite-size calculations. 



3 FERMIONIC REPRESENTATION 

At r = 27 (7 = 1, 2, ... a positive integer), the 2D jellium with the interaction 
Boltzmann factor rij<fc \ z j ~ z k\ F 1S mappable onto a discrete ID fermionic 
theory |10| . The mapping can be readily extended to the present model. The 
partition function l|2.7|l is expressed as an integral over two sets of Grassman 
variables V'j Q ''}j each with 7 components (a = 1, ...,7) defined on a 

discrete chain of N sites j = 0, 1, . . . , N — 1 as follows 



VtpV^e s ^ (3.1) 



7 (JV-1) 

S(£,l/0 = E E i w ^k (3.2) 

j,fc=0 

Here Vi^V^ = Y\^=q dip^ . . . dV'^d^ 7 '' . . . d£^ and the action S involves pair 
interactions of "composite" variables 

N-l N-l 

Bi= E 4 1) -i 7) - ** = E (3-3) 

Jl 1 • • • 1 j*y = fc]^,...,fc-y=0 

Ui+—+i~f=i) (*!+...+%=*) 
The interaction strengths w^fc [j, fc = 0, 1, . . . , 7(iV — 1)] are given by 

Wjk= j d 2 zw(z,z)exp[^jzjexp(^kz] (3.4) 



A 



Using the notation (• • •) = J 2>02?£e s • • • /Zn for averaging over the anticom- 
muting variables, the particle density (|2.1Q(I is expressible in the fermionic form 
as follows 

n(z,z)=w(z,z) ( E 3^k)eyLv(—jzjexp(—kzj (3.5) 



The fermionic correlators {(Sj^k)} can be obtained from the partition function 
<|3.1[l using relation l|3.2|l via the derivatives 

(S i * Jb ) = -^-lnZ JV (3.6) 
ow jk 
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Inserting the one-body Botzmann factor w of interest (|2.8[) into H3.4[) . and 
using the orthogonality relation 



f w/2 f 2vr 1 

J-W/2 ^ e ^\W i{j ~ k)y \= W6 > k 



(3.7) 



rW/2 
-W/ 

leads to the diagonalization of the interaction matrix 

1 f L/2 ( 2ir 1 

lOjfc = WjSjk, Wj = — j dx exp I -2-Kjnx 2 + — [2j - j(N - 1)] x > 

Notice that the diagonal interaction strengths {wj} possess the symmetry 

Wj = w 7(W -i)-j for all j = 0, 1, ... , 7(iV - 1) (3.9) 

With Wjk of the form (|3.8|l . the action 1|3.2[1 of the partition function (|3.1fl 
becomes "diagonalized" , 

A(iv-i) \ 
Z w = / P^X>£exp ^ E, (3.10) 

and the fermionic correlators {(Sj^fc)} are given by 

(H^ fe ) = (H;^)^, (H,* 3 -) = ^- lnZjy (3.11) 
The particle density (|3.5|l takes the form 

n W = W2 E ( S J*i) ex P -2^ 7 nx 2 + — [2j- 7 (iV-l)]x[ (3.12) 

where -AiV/2 < x < AiV/2, A being defined by Eqs. JOJ and lf2~Hjl . The 
symmetry (13.91 of the interaction strengths {wj} implies an analogous symmetry 
for the fermionic correlators 

(3,-*,-) = (H 7(JV _ 1) _ j * 7( ^_ 1) _ J .) for aU j = 0,l,..., 7 (AT-l) (3.13) 

This relation ensures the mentioned reflection property of the particle density 
n(x) = n(— x). 



4 GENERAL ANALYSIS 

The explicit formula l|3.12|l for the density profile contains the unknown set of 
fermionic correlators {(Sj\l/j)}. In this section we present an analytic treatment 
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of the thermodynamic N, L — > oo limit of this formula, supplemented by the 
periodicity assumptions l|2.12() and (|2. 13|) for the density functions with the 
period A defined by Eq. (|2.14|) . It will be shown that for any 7 the periodicity 
assumptions determine uniquely the general Gaussian structure of the density 
profile, without having at one's disposal the particular values of the fermionic 
correlators {(Ej^j)}. The treatment depends technically on whether 7 is an 
even or odd integer. 



4.1 7 = even integer 

For 7 even, we define the auxiliary integer M as follows 

7(iV-l)=2M (4.1) 

Let us shift the j-index enumeration in Eq. 1)3. 12[) by M, 

j = M + l, l = -M,-M+l,...,M (4.2) 

Now, rescaling appropriately the integration a;- variable in Eq. (|3.8|1 . the inter- 
action strengths can be written as 

exp(^) 1 f V%» f / /2F\ 2 1 

wm+i = 7== t= I ax exp < - \ x - W — / > (4.3) 



Here, we have introduced the dimensionless parameter 

(j, = nW 2 (4.4) 

which measures the number of particles in a square of side W; the limits fi —* 
and fx — > 00 correspond to the extreme ID (at zero temperature) and 2D versions 
of the model, respectively. The symmetry (|3.9|l is equivalent to 

wm+i = wm-i for all I = -M, ...,M (4.5) 

which can be verified directly from the explicit representation (|4.3|) . For / finite 
and in the limit N — > 00 

1 {2irl 2 \ , . 

wm+i ~ y= — exp (4.6) 



/2t7J V 

Under the index shift (|4.2|l . the density profile l|3.12|l can be expressed as 

— = J^- J2 c;(^)exp{-— (t--) 2 \ (4.7) 
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where we have introduced the coefficients 



ci(M) = (s M +i^M+i)J 7 r ex P\ — ( 4 - 8 ) 

The symmetry of the fcrmionic correlators l|3.13fl is equivalent to (Bm+i ^m+i) = 
M-i) ■ The consequent relation 

ci (M) = c-i (M) for all I = -M, M (4.9) 

ensures the reflection symmetry n(x) — n(—x). We see that, even for a finite- 
size system, the particle density l|4.7ll results as a superposition of Gaussians, 
localized equidistantly at positions x = \l/~f (I = —M, . . . ,M), with the uni- 
form variance a 1 = X 2 /j, / (Attj) but with different position-dependent ampli- 
tudes. The value-structure of these amplitudes simplifies substantially in the 
thermodynamic limit discussed below. 

All formal algebra made till now was rigorous. To describe the thermody- 
namic limit of the profile relation (|4.7|l . we adopt the two assumptions presented 
at the end of Section 2. In the limit M — > oo, one has to distinguish between 
the subsequences of M in (|4.1I) which correspond to N =cvcn and to N =odd 
particle numbers. Within each of the even and odd M-subsequences, the coeffi- 
cients {ci(M)} tend uniformly to their asymptotic values denoted by {c| e ^} and 
{c[°' 1 }, respectively. Thence, 

The analogue of the symmetry relation (|4.9|) takes the form 

c (e,o) = c (e,o) fo ra in = 0,±l,... (4.11) 

The periodicity assumption (|2.12() implies 

c (e,o) = c (e,o) f or a ii Z = 0,±1,... (4.12) 

The shift condition l|2.13|l between the "even" and "odd" states leads to the 
relations 

c\ e) = for all I = 0, ±1, . . . (4.13) 

Based on Eqs. (|4.11|) - (|4.13|) we conclude that there exist 7/2 + 1 independent 
asymptotic amplitudes Co, C\, . . . , Ct such that 

4 e) =C h c, (o) =Cj_, for 1 = 0,1,..., I (4.14) 
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All other coefficients can be generated from the basic set (I4.14|) by using the sym- 
metry relations 14. lip and 14.12(1 . The values of the asymptotic C-amplitudes 
depend on the dimensionless parameter fi of Eq. (|4.4|l . They are constrained by 
the neutrality condition l|2.11[l , written for x ranging over one period as follows 



d.r 







Simple algebra gives 



C +2{C 1 + --- + C, 1 ) + Cx=l 



(4.15) 



(4.16) 



For instance, in the 7 = 2 case, there are two amplitudes Co and C\ such 



that 



» 



Co for / even 

Ci for / odd 

Cl for I even 

Co for I odd 



The amplitudes are constrained by 

CoOO + dOiO 



(4.17) 
(4.18) 

(4.19) 



4.2 7 = odd integer 

For 7 odd, the fermionic formalism needs to be developed separately for the 
thermodynamic-limit subsequence with TV even and N odd. 
When the particle number N is even, the product 

7 (N - 1) = 2M + 1 (4.20) 

with some integer M. The shift of the j-index enumeration in Eq. I|3.12|) 



3 = M + - + 1, 



leads to the interaction strengths 



M 



-M ■ 



,M ■ 



exp 



w M+i+i 



(*) 



and to the coefficients 



N 



dx exp 



7M 



ci(M) = (E M+ i +l V M+ i +l ). — exp 

2 2 V 2[i \ 7/i 



(4.21) 



(4.22) 



(4.23) 



10 



possessing the symmetry q(M) = c_;(M). In the limit N — > oo (keeping AT 
even) , one arrives at the density profile of the "even" state 

The asymptotic coefficients satisfy the symmetry relation 

c\ e) = c^] for all I = ±i, ±|, . . . (4.25) 
and the periodicity condition 

c ; (e) = cg 7 for all I = ±i ±|, . . . (4.26) 

implied by Eqs. (|2.12|l and l)2.14[l . As a consequence, there exist (7 + l)/2 
independent amplitudes Ci , C|, . . . , Ca such that 

Ci (e) =Q for J = ±,§,...,2 (4.27) 

All other coefficients can be generated from the basic set (I4.27|) with the aid of 
the symmetry relations H4.25JI and l|4.26|l . 

When the particle number A" is odd, the product j(N — 1) is expressible 
similarly as in Eq. Q4.1JI . so that we can proceed along the lines of Section 4.1. 
The density profile of the "odd" state is again of the form (|4.10|) 

where the asymptotic coefficients exhibit the symmetries c, — c^°J and cf^ — 



c. 



1 — u -i 

;±7 



for all I = 0,±1,.... The shift condition (|2.13|) . when applied to the 
representations (|4.24|) and l|4.28|l , gives cf^ = cw_2. ■ In view of Eq. I|4.27|l , the 

basic set of the odd coefficients is given by 

4 o) =Cx^i for 1 = 0,1,...,^ (4.29) 
The neutrality conditions of type (|4.15(l constraint the C-amplitudes as follows 

2[Ci(m) + ---+Cj_ 1 (/*)]+C3(aO=7 (4.30) 
For 7 = 1, one has the simple result 



c 



l x xui oj± c _i_ 2 , _i_ 2 

» 



1 for alU = ±i,±|,... (4.31) 
1 for all/ = 0,±1,... (4.32) 
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For 7 = 3, the above scheme results in 

c (e) 

(o) 
c l 

where k is an arbitrary integer. The ^-dependent amplitudes are constrained 

by 

2Ci(/i)+C , |(/x) =3 (4.35) 

5 FINITE-SIZE CALCULATIONS 

In this section, we check the obtained analytic results by the exact finite-iV 
calculations. The crucial problem is to determine the dependence of the parti- 
tion function Zff{pf), given as the integral over anticommuting variables by Eq. 
<|3.10[1 . on the set of interaction strengths . Having at one's disposal 

this dependence, the consequent fermionic correlators [generated by using Eq. 
<|3.11[l ] determine the c-coefficients of interest via relations (14 .8f) [(7 even, N 
arbitrary) or (7 odd, N odd)] and I4.23|) (7 odd, N even). 
For 7=1, one has the simple result JT] 

Zjv(1) = w wi ■ ■ ■ w N -i (5.1) 
= — , 3 = 0,1,..., N -1 (5.2) 

Wj 

In the thermodynamic N — » 00 limit, after simple algebra both even l|4.31[) and 
odd (|4.32l) types of the c-coefficients are reproduced exactly. 

For larger values of 7, the partition function is a more complicated function 
of interaction strengths whose complexity increases with increasing the particle 
number N. The methods for a systematic generation of Zn(^), realized in 
practice through computer languages like Fortran, are summarized and further 
developed in Ref. ^2 • With the aid of these methods we were able to go up to 
N = 12 particles for 7 = 2 and up to N — 9 particles for 7 = 3. 

5.1 Results for 7 = 2 

The finite-size results for 7 = 2 are summarized in Figs. 1-3. Their discussion 
follows the analysis of Subsection 4.1. The asymptotic amplitudes Co and C\ 
are defined by relations (|4.17jl and H4.18JI . 

Taking the parameter [i = 3, the density profiles for the subsequence with 
the particle number N even and the one with TV odd are plotted in Figs, la and 



for I = ±A + 3fc 



for / 



3fc 



for / = 3k 

for Z = ±1 + 3* 



(4.33) 
(4.34) 
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lb, respectively. In the units of the period A = 1, hard walls are localized at 
x = ±N/2 where the boundary effects dominate. On the other hand, close to the 
center x = 0, the particle density exhibits the characteristic periodic behavior as 
from relatively small particle numbers. The density at x = has the minimum 
for N even and the maximum for N odd. This means that the asymptotic 
amplitudes fulfill the inequality Co < C\. The Gaussians with the smaller 
asymptotic amplitude Co are localized at the minimum points of the density 
profile, and so they modify the density plot given by the equidistant array of the 
Gaussians corresponding to the larger asymptotic amplitude G\ (and localized at 
the maximum points of the density profile) only marginally, without generating 
new extreme points. We observe numerically this phenomenon for any value of 
M- 

Taking the same value of the parameter /i = 3, the plots of the coefficients 
{co, ci, C2, C3}, reflecting the amplitudes of the Gaussians close to the center, 
versus 1/N are pictured in Fig. 2 where only the subsequence with N odd 
is considered. As the particle number N increases, cq — > C2 and c± — > C3 
as is expected. As can be seen the convergence of the c-coefficients to their 
asymptotic Co and C\ values turns out to be fast. 

The dependence of the asymptotic amplitudes Cq and C\, constrained by 
Co + C\ = 2, on the parameter fx is shown in Fig. 3. The finite-size errors in 
determining these amplitudes are deducible directly from the differences c — c 2 
and C\ — C3 at the largest N = 12 particle number. The error bars increase with 
increasing /i, they are however still so small that we do not present them in the 
figure in order to keep the clarity of the presentation. In the limit (j, — ► 0, which 
corresponds to the ID jcllium at zero temperature, Co = and C\ — 2; there 
exists consequently only one set of Gaussians (more precisely, the <5-function 
peaks), in agreement with the exact ID solution [3]. In the limit fi — > 00, which 
corresponds to the bulk 2D jellium, both asymptotic amplitudes tend to unity 
as it should be. 

5.2 Results for 7 = 3 

The finite-size results for 7 = 3 are summarized in Figs. 4-6. We do not 
comment these figures since they provide similar information as the previous 
ones for 7 = 2. We only recall that the asymptotic amplitudes C1/2 and C3/2, 
constrained by 2C\/2 + C3/2 = 3, are defined by Eqs. (|4.33|) and 14.3411 . 

One concludes that the finite-size calculations confirm with a high accuracy 
the predicted Gaussian structure of the crystalline state. The structure of the 
crystalline state for noninteger values of 7 is an open problem. 
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Fig. 6. The plots of the asymptotic C-amplitudes vs. \x: 7 = 
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